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Abstract. We consider electro-optical constructions in which the Casimir force is 
modulated in opposition to piezo-crystal elasticity, as in a stack of alternating tunably 
conductive and piezo layers. Adjacent tunably conducting layers tuned to conduct, 
attract by the Casimir force compressing the intermediate piezo, but when subsequently 
detuned to insulate, sandwiched piezo layers expand elastically to restore their original 
dimension. In each cycle some electrical energy is made available from the quantum 
zero point (zp). We estimate that the maximum power that could be derived at 
semiconductor THz modulation rates is megawatts/cm"^ ! 

Similarly a permittivity wave generated by a THz acoustic wave in a single crystal 
by the acousto-optic effect produces multiple coherent Casimir wave mode overtones 
and a bulk mode. We model the Casimir effect in a sinusoidally graded medium finding 
it to be very enhanced over what is found in a multilayer stack for the equivalent 
permittivity contrast, and more slowly decreasing with scale, going as the wavelength 
Acoustic waves give comparable theoretical power levels of MW/cm"^ below 
normal crystal damage thresholds. Piezo thermodynamic relations give conditions 
for effective coupling of the Casimir bulk mode to an external electrical load. Casimir 
wave modes may exchange energy with the main acoustic wave too, which may partially 
account for THz attenuation seen in materials. We outline feasibility issues for building 
a practical crystal power generator. 



PACS numbers: 77.65.Dq, 78.20.hb, 42.50.Pq, 43.25.Ed 
1. Introduction 



Though it is well accepted conceptually that the lowest energy quantum state contains 
energy, the quantum zero-point (zp) or vacuum fluctuation energy seems regarded 
somewhat distrustfully by many physicists as a real energy. The zp energy seems 
indispensable in the self-consistent formulation of many different physical processes: pair 
creation and annihilation, spontaneous emission, the Lamb shift, electron fluctuations, 
black-body radiation, van der Waals particle, and vacuum cavity forces (see Milonni's 
excellent overview book [1] ) , but the possibility that it could provide a significant energy 
transfer for human use gets too close to looking like a violation of energy conservation, 
getting something for nothing, or perpetual motion. Though zp energy is well accepted 
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as the source for the synthesis of the universe out of nothing in the big bang by pair 
creation [2] with enhancement effects possible in dynamical expansion processes [3], its 
importance in more common physical processes seems less acceptable. Still methods 
that may give a significant transfer of energy out of the quantum vacuum for practical 
utilization have been shown to be theoretically possible [1] , and possible implementations 
are being discussed and explored [3 [6l [7]. 

Vacuum cavity forces, especially in the example of the Casimir effect, give a 
macroscopic manifestation of the zp energy. Casimir [HI [9] showed that large-area 
closely spaced parallel conducting plates in a vacuum exhibit a force of attraction, 
which arises due to the distribution of electromagnetic fluctuations in the quantum 
zp field introduced by the conductors. We might suppose that the Casimir force is 
conservative like gravity, and parallel conducting plates should lose in their separation 
whatever mechanical energy was obtained by their attraction. However, it seems that 
this is really not the case: We know that the Casimir force can be turned on or off 
as tunably conductive plates are changed to conduct or insulate. If parallel tunably 
conductive plates operate elastically as against a spring between them, the work done in 
the attraction of the plates when the Casimir effect is turned on need not be returned in 
their subsequent separation when the Casimir effect is turned off. That excess of energy 
might be harvested, in effect, taking energy out of the quantum zp. The modulation of 
the Casimir effect by the conductive cycling of semiconductors with laser light has been 
demonstrated in precise differential tests, which have validated the theoretical estimates 
to within 1% [101111]. 



electrode 




Figure 1. Thought Experiment: Multilayer stack of alternating tunably conductive 
and piezo layers. Changing the tunably conductive layers to conduct and then insulate 
produces a Casimir compression followed by its elastic reexpansion, which does net 
work on the piezo crystals from which useful electrical energy might be derived. 

To understand the possibilities for energy extraction, we consider as an idealization 
a stack of alternating tunably conductive and piezo layers bounded by conducting 
electrodes as illustrated in Figure [1] The Casimir effect between two conductive layers 
produces a net attractive force per unit area F in the intermediate piezo layer that 
drops off rapidly with the increasing piezo layer thickness £-, so the piezo layers must 
be quite thin to obtain a strong effect. Sandwiching conductive plates separated by a 
small intermediate dielectric of permittivity Ee- and thickness i- exhibit a near-field 
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Casimir force F(^{i_) that follows al/i^ power law up to some cross-over separation i^, 
and then switches to a steeper power law, written for perfectly conducting infinite 
area bounding layers from the Lifshitz formulation p[2] 



Tc^hc 



270 nm 



1 4 



dyn/cm^, (1) 



where the static permittivity of the piezo layer usual h is Planck's constant 

over 271 and c the speed of light. For the rightmost equality, a static permittivity of 
ee-=6 (index-of- refraction squared) was used for moderate- index insulators, and gives a 
Casimir force per unit area of just 1 dyn/cm^ or 1 erg/cm^ = IE— 7 Joules/cm^ for the 
piezo thickness = 270 nm (using the 'E' notation for the power of 10). Though the 
Lifshitz solution is for infinitely thick sandwiching plates, his solution represents a good 
approximation too when the plates are comparable in thickness to their separation 
~ i_, as we discuss. The cross-over separation where the power law switches to 
is characteristically less than 7 nm for near conducting sandwiching layers, as we 



elaborate in Appendix A 



In acting against an intermediate elastic material, the attractive force per unit area 
F produces a dimensionless strain —sF or a compressional displacement sFi_ in a layer 
of thickness where s is the material elastic compliance for the stress component in 
the forcing direction. The average work done per unit area by the Casimir force in 
compressing the intermediate layer by a continuous motion is the integral through a 
half cycle of the force F times the infinitesimal displacement si^dF or si-F'^/2 for F 
the maximum force in the cycle, giving the total work done per unit of time and volume 

W = ^suF\ (2) 

for cycling at a rate u, and supposing equal odd-even layer thicknesses = so 
half of the volume participates in the compression at any time. This formula is an 
approximate form useful for estimating the amplitude of the effect, as it supposes the 
idealization of high piezo efficiency, a static piezo response, and ignores contributions 
from the converse piezo effect, as we elaborate in Section [61 

Semiconductor modulation at rates up to about z/=lE13 Hz by electrical 
modulation in adjacent oppositely doped layers for modern transistor design has been 
demonstrated [13]. Material compliances range from about s ~1E— 7 to IE— 12 Pa~^ for 
soft rubber to the hardest crystals like diamond. For estimating the amplitude of the 
effect, we suppose most conservatively s ~ IE— 12 Pa~^ = IE— 13 cm^/dyn, and obtain 
a power per unit volume 1^=0.25 erg/s/cm^ in a stack with layer thickness i- = 270 nm 
from (|2]). As the power scales inversely as i^, with a tenth that layer thickness £_=27 
nm, we obtain a much larger W=2.5 Watts/cm^, and for ^_=5 nm, W=1.8 megawatts 
(MW) / cm^ ! Actually the estimate is somewhat high for normal good conductors, which 
fall significantly short of the ideal of perfect conductivity. Still our more careful estimate 
gives power up to megawatts/cm^ for idealized stacks of alternating tunably conducting 
and piezo layers, as we discuss in Section HI Physical constructs could be possible as we 
elaborate in the text as deposition of layers only £_=5 nm, about 10 molecular layers 
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^ ^ 




Figure 2. ZP energy from a standing acoustic wave in a piezo crystal. Successive 
close permittivity plane wave fronts introduced by an acoustic wave interact like 
alternating layers, giving an added Casimir force that is modulated in time with the 
wave oscillation. 

From the thought experiment, we recognize a basic physical principle and study 
processes that systematically modulate the Casimir effect, like a pressure wave standing 
or traveling, longitudinal or transverse, which produces a periodic sinusoidal permittivity 
variation by the acousto-optic effect, as illustrated in Figure O As with the alternating- 
layer stack, adjacent plane wave fronts in permittivity must interact via the Casimir 
effect giving an added longitudinal stress at the extremes of the temporal oscillation, as 
illustrated by the force vectors F. The Casimir effect in acoustic waves too may produce 
a voltage difference across a piezo crystal from which electrical power might be derived. 

The Casimir effect in graded permittivity media has been considered mainly to 
understand effects of soft boundaries on three-layer solutions, and shown to exhibit 
quantum divergences unlike what is found for solutions with discrete layers, suggesting 
an enhanced effect O HSl |16]. We formulate the problem again for a general ID 
spatial graded permittivity variation in Sections [2] and El Numerical models for stacks 
of alternating permittivity layers are presented in Section HJ and for permittivity waves 
in uniform media in Section |5l Besides being enhanced for the equivalent permittivity 
contrast, the Casimir effect also exhibits a much slower falloff with the scale of the 
variation than what is seen with discrete layers. 

The zp energy dynamics are examined using the piezo thermodynamic relations in 
Section |6l which shows formally the transfer of energy out of or back into the quantum 
zp in the Casimir modulation of piezo crystals for certain acoustic wave conditions. 
Energy that is not used is cycled into and out of the quantum zp, with the strain and 
stress modulated in phase. Estimates for power production for alternating-layer stacks 
and acoustic waves in crystals from our numerical models are given in Section [71 Casimir 
wave modes may interact with the principle acoustic wave too to drive its attenuation 
or growth, as we discuss in Section |H1 Some practical issues for design of a zp-power 
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harvester based upon acoustic waves in single uniform piezo crystals are discussed with 
conclusions in Section |9l 

2. Casimir Effect with a ID Permittivity Variation 

The Casimir force arises in a medium with a ID permittivity variation z) in a spatial 
coordinate z as in the multilayer stack illustrated in Figure [1] or the single crystal in 
Figure El The general complex permittivity z) characterizes fully the electro-optic 
properties of a medium, defining the index-of-refraction and absorption variations as a 
function of light frequency ^. The Casimir force is derived in many places especially 
after the general Lifshitz formulation [12]. We write it as a force per unit area or a 
pressure as a function of z with positive as compressive 

F{z) = ^ V ' / tuK{z) {R,{z) + Rp{z)) dzu, (3) 

for /cb Boltzmann's constant and G the temperature. The sum is over the Matsubara 
frequencies ^ = n27r/cB0/^ = r;,[6/300K] 2.47E14 rad/s, where the integer index n 
counts the frequency intervals, and the prime on the summation indicates that the 
zero frequency n = term receives half weight. The two terms R^a^^ represent the 
Fresnel reflection variation for perpendicular s (senkrecht) and parallel p zp polarized 
rays deflned 

R^s^^+{z)R^sy{z) 

for Fresnel reflection coefficients which depend upon the permittivity. The 

wavenumber n{z) is for zp dissipative electromagnetic modes and is deflned in by the 
Matsubara sum and integration frequencies ^ and w 

n{z) = -^{e{i,z)e + w^y^ . (5) 
The deflnition for k.{z) and Fresnel reflection coefficients R(^s^j^ use the function 

e{E,i z) = i{iC,, z), the projection of the complex permittivity on its imaginary frequency 
axis, which casts the complex e into a real function e consistent with Kramers- Kronig 
causality constraints [E]. This solution may be applied to represent index, absorption, 
conductivity, or magnetic permeability variations in a medium. 

It is usual to replace the finite sum over Matsubara frequencies in (j3]) by an 
integral. Symbolically the sum might be approximated Yl'^=o ' ■ ■ ■ ^ /o°° • • • dn, which 
takes properly into account the half-size interval around the lower limit = 0. Using 
the Matsubara formula with = {2TikBQ/h)dn gives on substitution the integral 
approximation for the Casimir force 

^TT C Jq Jq 

A general substitution for the sum that takes into account singularities in the complex 
permittivity function i{i^) must include Cauchy principle values as represented by the 
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Abel-Plana formula [IB]. The singularities in the complex function e{i^) are on the 
positive real axis in its complex argument iC, and enclosed in the upper half plane of i^ 
by the infinite integral along the imaginary axis of i^, as contained in the sum over the 
real Matsubara frequencies ^ from zero to infinity in the formula ([3]). Simply substituting 
the integral for the sum represents a good approximation for all temperatures when the 
scales of variation are smaller than about 1 /xm [I9] as we illustrate in numerical models 
in Section m 

Fresnel reflection coefficients R^s'^^{z) for the relevant dissipative waves of 
wavenumber k{z) in the s and p polarizations are propagated by recursion on the left 
side (— ) from R^.s-^_{z) = as 2 — )■ — oo or on the right side (+) from R^s-^_^{z) = as 
z — )► +00. The left-side reflection coefficients (o^ right-side R^=^_^_) change moving 

a distance Az (or moving —Az for the right-sided coefficients) to cross a permittivity 
jump at some z following the comprehensive formula [Section 1.6] [20] 



r(s^Jz) + Rfssjz ± Az)e-2«(^±^^)^^ 

p ( ^] = VP/ /y\ 

l + r^s^^{z)R^s^^{z±Az)e-M-±^z)Az^ y' ) 

with Az always positive, where the coefficients for the jump at z, r^^s^^^z), are defined 
by the relative differences 

k{z =1= 6z) — k{z ± 6z) . . 

^ K{zTSz) + K{z±Szy ^ ^ 

e{z ± 6z)k{z ^ 6z) — e{z ^ Sz)k{z ± Sz) . , 

""^^^""^ " e{z ± 6z)ti{z T Sz) + e{z t Sz)k{z ± 6z) ' ^ ^ 

for any infinitesimal (positive) distance 6z < Az. We observe the sign convention that 
the differences are between the wavenumber nearer z = minus the wavenumber further 
out, which leads to convenient solution symmetries. With no permittivity change over 
an interval Az, the jump coefficients are zero f^^"j± = 0, and the Fresnel reflection 
coefficients decrease exponentially in +z for the left-sided solutions R^s-j_{z) or in —z 
for the right-sided solutions i?^s^^(z) with i?^s'j_,_(z) = R^s^^{z ± Az) exp{—2KAz) . The 
formulae are relevant for reflection coefficients crossing discrete permittivity jumps, but 
are applicable in graded-media too, as we elaborate in the next section. We sometimes 
denote the variables as £^(0 or Kj using subscripts rather than functions of z for discrete 
layers counted from a middle layer j = as illustrated in Figure [3l The ± subscript 
is also used for even and odd layers with alternating-layer stacks, where the reversed 
state is of special interest with interchanged permittivities ^ and wavenumbers 

K,-^ Kj— . 

For a symmetric permittivity variation e{z) = e{—z), recursive application of 
the formula ([7j) gives R^s^^j^{z) = R^s-^_{—z), or for the product R^s-j_^{z)R^s^_{z) = 
Rl^s'^_^{—z)Ri^s^_{—z), for Jli^s^^{z) = il^s^(— 2;) from (jlj), and a symmetric Casimir force 
in © or (ED with F{z) = f(-z). 

In a multielement stack or square permittivity wave, the Casimir force can be seen 
to be constant in every constant permittivity zone. The left /right-sided (=f) Fresnel 
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Figure 3. Plane-parallel multilayer sandwich with layers of differing permittivities 
on each side of a layer j = of principle interest 



reflection coefficients are propagated from bounding permittivity jumps at z = =F^o/2 
to a given 2; in a middle layer of thickness £0 as in Figure |3] with 

i?(.)^(^) = i?(.)^(±|)e-^-(^^^), (10) 

so in the product, the z dependence vanishes with 

= fl(.),(|)fl,.)J-|).-'., (U) 

giving H^s-j from (jlj) and F from ([3]) and ([6]) independent of z across the layer. 

Solution z independence is evident across the middle layer in the multilayer 
solutions given in the literature. For a three-layer stack the jumps at 
z = ±£0/2 define the products of Fresnel reflection coefficients R^s^^R^s^_ = 
r^3^^(£o/2)r^s^_(— £0/2) exp(— 2ko£o) from ([8]) and Qj, which agrees with the Lifshitz 
formula (Eq. (4.14) in Dzyaloshinskii et al. |21j, with some substitutions). For 
a five-layer stack = 2 with layers of thicknesses ^j counted by —2 < j < 2 
as in Figure [HI we obtain products of Fresnel reflection coefficients as in ( fTTj) with 
the reflection coefficients -Rj^s-j_,_(±£o/2) defined by (171) using the two upstream jump 
coefficients r^s^_,_(±(£-|-i + i'o/2)) defined in ([8]) and ([9]), which agrees with Zhou and 
Spruch (Eqs. (3.8) and (3.12)-(3.16) in [22]). The formulae give similar agreement for 
the Casimir force in ([3]) or ([6]) with any number of layers as developed in numerous 
studies [231I21I2S]. 



3. Casimir Effect in a Graded-Permittivity Medium 

The definition for the Casimir force from Eqs. ([3]) or ([B]) is applicable too in a 
graded-permittivity medium with a ID continuous permittivity s{z), as can be seen 
by approximating the permittivity variation by discrete layers, for some small layer 
thickness Az, as illustrated in Figure |H A spatially shifted form of the recursion 
formula ([7]) lends itself directly to numerical computations for propagation of the Fresnel 
reflection coefficients far from the domain of interest with initial = through 

intervals with centered permittivity jumps written for the left-sided coefficients 
r/.N (2 - Az/2) + R(s^{z - Az)e-'^(^-^^)^^ 
(;)-^ ^ l + r^sy{z-Az/2)R^sy{z-Az)e-<'-^')^' ' ^ ' 
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for jump coefficients ri^s^_{z — Az/2) still defined as in (jS]) and (jH])- The right- 
sided reflection coefficients could be propagated similarly from z — )■ oo in 
the —z direction, but in our numerical models they are defined using the relation 
+ = applicable to symmetric arrangements around a middle z = 0. 
A direct computational approach like this one has been applied to estimate the Casimir 
effect in linearly graded media by Inui [26] . 




Figure 4. Direct calculation of Fresnel reflection coefficients in a graded-permittivity 
medium. A continuous permittivity variation e{z) is approximated by discrete layers 
of small thickness Az. 



The Fresnel recursion formula for a graded-index medium ([7]) or ( !T2|) with (|8]) and 
dH]) is written in the limit that the interval Az becomes infinitesimal Az — ?■ by the 
differential equation for the continuous variation of the Fresnel reflection coefficients 
R^s^^{z), as in the example [SOl Section 1.6]. Retaining only terms to ffist order in Az, 
we obtain from either recursion relation 

with the jump coefficients approximated by derivatives for the spatially continuously 
permittivity e(z) 



2k{z)Ris\_{z^, 



(13) 



2k(z) dz 
1 dn 



2ck[z 
1 de 



de 
dz 



2k(z) dz 2e(z) dz 



2ck(z) 




(14) 
(15) 



with the wavenumber k{z) still defined by Eq. ([5]); e{z) and k,{z) retain their 
dependencies on the frequencies ^ and/or zu implicitly. 

We examine bounded solutions with everywhere- small reflection coefficients 
l-R(^s^_(-2)| ^ 1; which is consistent with our main interest in small periodic permittivity 
fluctuations on top of a constant permittivity background. Thus we are able to linearize 
Eq. f|T3l) and obtain the closed-form integral solution 



R^syiz) 



3— 2 f K(z)dz 



^2f K.{z)dz 



1 - i?/, 



dz 



(16) 



The function R'(s\ (z) represents a prior approximation to the Fresnel reflection 
coefficients, which we apply in our numerical iterations. For our analytic evaluation. 
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we leave the integrals all as indefinite and functions of z. A lower bound on the 
integral in the exponential integrating factor is of course arbitrary, as any constant 
in the exponential factors out of the main integral to cancel with the same constant 
from the outer exponential. 

With an added constant of integration in the square brackets in ( IT6l) . a term is 
added to the reflection coefficient due to the leading external exponential integrating 
factor, which diverges as z — )■ — oo for Rf^s^^^^z) as k{z) is always positive, or as z — )■ +00 
in the similar expression for R^s-^j^{z). Leaving the main integral in square brackets as 
indefinite always gives bounded solutions. 

Bounded refiection coefficients must also be periodic with a periodic permittivity 
fluctuation, as shifting the origin of the calculation by any number of wavelengths nX 
in z cannot change the form of the solution. For a periodic permittivity fluctuation, the 
wavenumber from (|5]) can always be written as a constant plus a spatially periodic 
term k{z) = + lii{z), giving an exponential integrating factor that decomposes 
into the product of constant and spatially periodic functions exp(2 f n{z)dz) = 
exp(2Ko2;) exp(2 J ki{z)dz), as the integral of a periodic function with zero average like 
ki{z) is also periodic. Thus shifting the origin of the main indeflnite integral in (fT6l) by 
an integer number n of wavelengths z ^ z + nX, factors the constant exp(2nft;oA) out 
from the square brackets leaving the square brackets unchanged, and this constant just 
cancels with its inverse from the outer exponential. Even the perturbation term R'/s\ (-2) 

deflned by a previous iteration must be periodic, so periodic solutions generalize beyond 
the linear case with R^s\^_{z) = Ri^s^^_{z + nX). 

Leaving the main integral as indeflnite in ( flGll is convenient for our analytic 
evaluation, but in numerical implementations, the integral must be evaluated between 
specifled limits. Then we include an added constant of integration in the square brackets, 
which we deflne by the periodicity condition R^s^_{z) = R^s-^_(z + X) for arbitrary z for 
periodic permittivity variations. A computationally desirable form for the constant of 



integration for numerical applications is derived in Appendix B 



The Fourier decomposition for the spatial variation of a permittivity wave in z with 
minimum permittivity at z = is written 

£:(z) = £0 - e^iccos— , (17) 

using the tilde ek{cs) to denote the (real) spatial Fourier wave mode with wavenumber 
k in cosine (c) or sine (s) elements. The Fourier wave modes retain the Matsubara 
frequency ^ dependence implicitly Eq = £o(0 cind iic = £ic{0- We do not introduce 
temporal varying factors here as wave changes are always assumed to be slow compared 
to frequencies that are important for the Casimir effect. The permittivity variation 
e{z) might represent a standing wave like what is indicated in Figure |2l which does 
not change its spatial form but exhibits a sinusoidal temporal variation through an 
external temporal factor, as e{z) — i- e{z) cos27rz/i(: or in the usual analytically continued 
form with e{z) — )■ £:(z) exp(z27rz/t) and cos(27rz/A) — )■ exp{i2TTz/X) as we employ 
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in spatial/temporal Fourier analyses in Section O Solutions for traveling waves are 
obtained similarly with a ±z moving wavefront by the replacement z — t- Xut =1= z. 

Assuming a relatively small wave amplitude \eic\ ^ \eo\, the jump coefficients 
r^s^_{z) and the wavenumber k{z) are written to ffist order in iic using the z derivative 
of the permittivity from ( |T^ and ( |T5i) and the definition 

27TZ , , 

sin—, 18 



= sm- 



= — Kic cos — — , (19) 



where 



^Sls „ 9 , _o 5 '"pis '"sis ~ ) 

2c^Xk,q Xeq 

f^o = ^(^0^^ + ^'^Y^^, and = 

Note that to ffist order, kq is independent of iic, and all of the other coefficients are 
proportional to e^ic. 

Substituting into ( IT6II and retaining terms only to ffist-order in the permittivity 
fluctuation iic, leaves the analytic integral for the Fresnel reflection coefficients 

9- /■ 9- 271Z ^(^)isKosin^ - ^cos^ , , 

R,^. (z) = r.^. e-^^°' / e^-^"" sin — d;^ = ^ ^ ^ ^. 20 

The indefinite form of the integral gives a periodic solution as we expect. 

The Fresnel reflection coefficients Ri^s\^j^{z) make contributions to the Casimir force 
in ([3]) or ([6]) through their products R^s\^^{z)R^a^_{z). Remembering that R^s^^{z) = 
R^s^_{—z) for a symmetric permittivity variation around z = as in the wave from 
(IT7|) . we obtain 



Rr 



87r2 

defining the new variable x = ^qX/ti. 

The jump coefficient fpis from ( IT5l) is of order and f^is from ( IT4l) of order or 
smaller than the relative permittivity fluctuation Sic/eq, being possibly smaller since 
Kg = eoC^/c^ + tu^/c^ > eq^'^/c^. Thus the product i?j^s-j_|_(z)i?j^s-j_(z) is no larger than 
second order in the relative permittivity change. The offset in the denominator of 
R^s^{z) in (jlj) introduces a still smaller-order effect, so for a small-amplitude standing 
permittivity wave, the Casimir force from Eq. ([6]) is written to most-significant order 
in the relative permittivity fluctuation as 

i^W = ^^ X X^Rs+{z)R,.{z)dxdi + x'Rp+{z)Rp~{z)dxda, (22) 

having eliminated the variables kq and by substituting with the new variable x- The 
jump coefficients from (iMl) and (fTSjl are written in terms of x as 

^ _ ^^Aeic - _ - _ ^£ic fr.o\ 
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The polarization product Rs+{z)Rs_{z) goes like l/x^ asymptotically for x 3> 1, 
which makes a contribution to the first integral in fl22l) . whereas the second 

polarization product Rp+{z)Rp_{z) goes like l/x^ at large X; so the second integral 
in ( 122|) is linearly divergent. The second divergent p polarization integral is dominant, 
so we drop the integral over s polarizations in (l22ll and obtain 

having factored out the ^ integral with the permittivity contrast Sic/sq, supposing 
no dependence upon ^ in the cutoff Xc- The x integral cutoff is determined by a 
cutoff frequency that is much higher than the maximum frequency for permittivity 
fluctuations, which is the electronic vibration frequency, and so can be approximated 
Xc ^ Awc/(vrc). 

The divergent integral is treated by introducing a smooth Lorentzian cutoff function 
into integrations around an upper Xc in our numerical calculations as a form of covariant 
regularization [271 Section 9.6]|28l Section 9.2], and by simply truncated the integral at 
Xc in this analytic derivation. Both procedures give very similar estimates for the 
Casimir force as we discuss in Section [5l 

Regularization by imposing a maximum frequency cutoff natural to the physical 
process seems sufficient for our interest in estimating the amplitude and in exploring 
general properties of the Casimir effect in permittivity waves. We argue that electronic 
reactivity to virtual photons must be limited by a maximum recoil due to the small but 
finite electron inertia, or perhaps equivalently the qed calculation exhibits a small but 
finite correlation length represented by the Compton wavelength [29l Chapters 8-9]. In 
either case, the Compton frequency seems the only relevant choice for an upper frequency 
cutoff in ordinary materials, as is used in other studies that test the Casimir effect, and 
as has been successfully applied with other zp effects, historically and most famously 
with the Lamb shift [30], [281 Section 9.6.2]. Thus we take ojc = rUeC^ / h=7 .7 6E20 rad/s 
for nie the electron mass. 

The X integrals in (12^ are analytic and both well approximated by the value of 
their upper limits Xc for Xc^ ^, thus giving 

The Casimir force in a medium with a small-amplitude relative permittivity 
fluctuation goes as at all wavelengths, and is second-order in the relative 

permittivity fluctuation. The force exhibits significant power in its spatially constant 
and cos(47r2/A) wave modes only, and ranges from zero, where cos = 1, to a maximum 
negative value where cos = — 1. Thus the Casimir effect in small-amplitude permittivity 
waves is always repulsive. We refer to the Casimir modes separately equating terms in 
the spatial Fourier expansion F{z) = Fq + F2cCos47rz/A. 

This largest-amplitude spatial wave mode contributing to the Casimir force from 
( l25l) does not change sign with the changing sign of the permittivity wave, as the 
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permittivity contrast enters only with its squared value. The Casimir force vanishes 
when the wave permittivity contrast vanishes at zero temporal phase, and so the Casimir 
force is modulated with twice the permittivity wave frequency 2z/. It can be seen that 
spatial wave modes other than cos(47r2;/A) exhibit terms proportional to ^(^^'^is and so 
do exhibit variations at the normal wave frequency z/, but these terms are of lesser 
order in powers of the relative permittivity fluctuation. We study the properties of the 
smaller-order Casimir wave modes numerically in Section |5l 



4. Casimir Effect in an Alternating-Layer Stack 



Figure [5] illustrates the Casimir force in the middle layer of an alternating-layer stack 
as a function of the spatial scale A = £_ + for permittivities e±(0 i^i the drawing 
Figure [1] or from Figure [3] with equal layer thicknesses i+ = The three-layer 

calculation {n^ = 1) derived by integrating numerically Eq. ([6]) is shown with three-layer 
calculations for three temperatures using the Matsubara sum from (j3]) {solid lines), with 
the 41-layer calculation (n^ = 20) with (I6l) {dot dot long-dash), and with limiting power- 
law forms {dot dash) and F) {long dash) from Eqs. flA.3|) and flA.6l) . respectively, 



discussed in Appendix A The lower panel illustrates the same data on an expanded 



scale with the 3-layer curve subtracted. 




Figure 5. Casimir force as a function of scale X — £- + £+, from IE— 8 cm=:lA to 
IE— 3.6 cm=2.5/im, with alternating static-permittivity contrast Asc/Sc = 0.01, for 
3 layers and for three temperatures (solid lines), for 41 layers Ui = 20 {dot dot long- 
dash), compared with limiting power laws F( {dot dash) and {long dash). The same 
data are shown in the lower panel on an expanded vertical scale with logj^Q F for 3 
layers subtracted. 



For our models, we adopt a usual form for the complex permittivity as a function 
of light frequency ^ 
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which is defined by a single electronic band at frequency oOq characteristically ~2E16 
rad/ s of strength 6^ — 1, and a conductive term defined by the plasma frequency Up, which 
ranges from to 2E15 rad/s in semiconductors, and reaches about 1.6E16 rad/s in doped 
semiconductors or in the best room-temperature metallic conductors, with the Drude 
collision frequency 7p, characteristically ~2E14 rad/s. The simplified form contains the 
important terms for the Casimir effect [17] . A lattice vibration band at lower frequency 
(typically 9E12 rad/s) represents the next lower frequency effect in crystals, but is both 
below the important frequency range and of relatively small oscillator strength to make 
a significant contribution to the Casimir effect. The main contribution to the static 
permittivity comes from the electronic vibration band with i{C, = 0) = Eq] ~6 for 
materials of interest. 

Due to causal restrictions, models for the complex permittivity must always be real 
on its imaginary frequency axis. The real projected permittivity from ( l26l) is defined 
e{^) = i{iC,), which gives 

The projected permittivity must be a slowly varying monotonically decreasing 
function [T7], and is plotted in many places for conductors or nonconductors (e.g. 
see Figure 20 [H]). For the calculations for Figure |5l the material is assumed to be 
nonconductive {up = 0), with a;e=2E16 rad/s and ec=6, with static permittivities Eej in 
successive layers j alternating between £e± = ± '^^e] for the figure the permittivity 
contrast used was A£:e/£^e=0.01. Following the convention used in earlier sections, the 
middle j = and even j layers are denoted by the subscript '— ' and the odd j layers 
by '+'. Though the formulation discussed is mostly applicable to arbitrary multilayer 
stacks with numbered layers, we remain interested primarily in stacks of alternating ± 
layers. 

For small permittivity changes |A£:| ^ e, the parametric representation in the 
model (l26l) and (1271) is assumed to remain accurate. Then the permittivity change can 
be described by the derivative expansion 

de , de , de , de , , , 

dee dbjp d-fp 

Physical processes in materials change permittivity parameters in different ways, so we 
explore the Casimir effect with different changing parameters in our numerical models. 
Although such parametric processes should give the most-significant contribution to the 
Casimir effect that can be introduced with acoustic waves, additional spectral effects not 
represented in the permittivity model of (126|) also enter. The permittivity spectrum e{C,) 
exhibits significant changes with pressure in the vicinity of lattice vibration resonances 
commonly described as a spectral variations in the photoelastic coefficient [311 Chapter 
9]. The photoelastic coefficient enters in the phenomenological acoustic wave theory, 
which we describe at the beginning of Section [HI We note in our numerical experiments 
the very different limiting forms that arise when different parameters controlling the 
permittivity are perturbed. 
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As is well-known and evident in Figure the Casimir force for a three-layer stack 
exhibits relatively small- and large-scale limiting power-law forms. Our Figure [5] in 
the 41-layer = 20 calculation {dot dot long-dash) suggests the generalization of this 
property to arbitrary multilayer stacks. As long as the jump coefficients for a multilayer 
stack are small Ir^s-j^,! ^ 1, their contribution to the reflection coefficients R(^s^^ in 
([7]) are increasingly exponentially attenuated with increasing distance from the middle 
z = 0. The final exponential factor in the recursion formula exp(— 2k(0)£o) as it appears 
in (fTTl) . must always be the dominant one, and increasing the number of layers cannot 
much change the limiting form of the integral. We reexamine and explore the scope of 



the Lifshitz limiting solution forms in Appendix A 



We find it convenient to use the spatial period in alternating layers as the relevant 
scale, i.e. the wavelength A = . Though the thickness of the middle layer £_ alone 

would be the natural choice, the full period facilitates comparison with the wavelength 
for the Casimir effect in graded media described in subsequent sections. Though we use 
wavelength for the dependent variable, and assume equal odd-even layer thicknesses in 
all the figures shown in this section with A = = 2i_, the Casimir effect in an 

alternating-layer stack is essentially determined by the even layer thickness i- alone, 
and changes little with changes in the relative thickness of the odd layers over 
several orders of magnitude. 

In our numerical calculations, the infinite-frequency integrals in ^ and w in Eqs. 
([3]) and are written as finite angular integrals by substituting for each ^ or zu with 
cjetan^^ and integrating from ^ = to 7r/2. The ^ integral is convergent because the 
permittivity goes to zero for ^ ^ Wg. Though the integral is divergent in general, 
an effect evident for graded media as we have discussed, it is convergent for multilayer 
stacks due to exponential attenuation of reflection coefficients from (ITT]) , which cuts off 
at wavenumbers k(0) ^ I/Eq in the middle layer or frequencies zu ^ k(0)c = c/Eq. 
We flnd that the integrals converge consistently in the fewest integration steps with the 
power L = Q, and obtain relative accuracies better than IE— 3 using Simpson's rule with 
120 steps, which is what we used for all the multilayer calculations described here. 

Figure |6] shows the Casimir force F{X) for alternating static permittivity Eg in (a), 
electronic vibration frequency cug in (b), and plasma frequency Up in (c). The static 
permittivity (a) and the electronic frequency (b) alternate with the small differences 
identifled in the flgure, but the plasma frequency (c) alternates between Wp = in 
the nonconducting even layers and Upj^ in conducting odd layers with 7p=2E14 rad/s. 
Results are overplotted for calculations made with increasing numbers of layers = 1, 
2, 3, 5, . . . (all shown as dashed), and = 20 {solid). Relative changes less than 5E— 3 
in the estimated Casimir force over the range of scales are found with increasing the 
number of layers used in the calculation above about = 3. 

The Casimir force estimates all follow the 1/A^ or 1/£q thin-limit power- law form, 
but the thick-hmit power law 1/£q does not apply with fluctuations of the electronic 
vibration frequency Aue/oOe as we elaborate in Appendix A The numerical experiments 
for the fluctuating electronic vibration frequency (b) suggest another deflnite thick- 
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Figure 6. Casimir force F{X = £_ + ell+ = in the central layer of a multilayer 

stack with alternating static permittivity Eg (a), alternating electronic vibration 
frequency cOg (b), and plasma frequency in the odd layers Wp+ (c). Calculations with 
different parameter contrasts for 41 layers ne = 20 (solid) and fewer n^=l, 2, 3, 5, 
. . . (dashed) are shown. 

layer power-law form of l/^g- Though power-law dependencies 1/£q and 1/£q are seen 
over ranges in A or with alternating plasma frequency in panel (c), these do not 
show a single cross-over separation, but effects depend strongly and specifically upon 
the plasma frequency contrast. 

The separation between the successive force curves with increasing permittivity 
contrast goes as the relative contrast squared, (Ase/^e)^ in panel (a), and as {AaJe/uJeY 
in (b), with deviations in the Casimir force across the range in A no larger than 5E— 3 
up to relative contrasts 0.01. For higher contrasts, larger variations arise as nonlinear 
effects begin to be important. A squared dependence is natural, as the Casimir force 
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from goes like . The reflection coefficients feel their main contribution from 

the first jump, and so must go like i?j^s-j_,_(£o/2) ~ '^(°)± discussion following 

(|TT1) . and the jump coefficients follow the relative derivative ^(^^"j± ~ {de/dz)/e in their 
most significant contributions to the integrals, as might be seen from (fT4l) and (ITSl) . 
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Figure 7. Casimir force difference AF(A) between normal and permittivity-reversed 
stacks for alternating static permittivity Ec (a), alternating electronic vibration 
frequency Wo (b), and plasma frequency in odd layers ujp+ (c). Calculations with 
different parameter contrasts for 41 layers ni = 20 (solid) and fewer ni=l, 2, 3, 5, 
. . . (dashed) are shown. 

Though the Casimir effect is mainly unchanged between the normal and the odd- 
even reversed stacks, small residual systematic differences do exist. The change in 
Casimir force between normal and reversed static-permittivity states AF{X) is shown 
in Figure [7] for alternating static permittivities in panel (a), alternating electronic 
vibration frequency Ue in (b), and plasma frequency in odd layers Up+ in (c). Force 
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differences have to be computed before frequency integrals in to avoid a large relative 
numerical noise, as the difference is small compared to the force itself. Force differences 
are overplotted for increasing numbers of adjacent layers used in the calculation = 1, 
2, 3, 5, . . . (all shown as dashed), and ne = 20 (solid). Relative changes no larger than 
about 0.05 are found for calculations made with more than about 7 layers {n£ = 3). 
Force differences between normal and reversed stacks AF{X) increase as the permittivity 
contrast cubed in (a) and (b), as is natural for small differences between forces that are 
the same to second order in their even and odd permittivity contrasts. 

Force differences exhibit good thick-layer power-law forms with fluctuating static 
permittivity Ee following the same power law as what is seen as for the Casimir force l/fg 
in (a), and a much steeper thick- layer power-law form l/^g" fluctuating electronic 
frequency Ue in (b). Force differences between a normal and reversed plasma-frequency 
change in (c) show characteristic power-law forms 1/£q and I/^q over ranges in scale, 
but with no obvious dependence of the cross-over separation upon the plasma frequency 
contrast. Force differences do not exhibit very well-defined thin- layer power- law forms, 
but go only approximately as l/^g^ in all panels (a)-(c). 

As the modulation of the Casimir force cycles between the nearly identical 
normal and reversed permittivity states it passes through zero, so with a permittivity 
modulation at frequency z/, the Casimir force is modulated mainly at the frequency 2i/. 

5. Casimir Effect in a Permittivity Wave 

Figure [H] shows an example Casimir force profile F{z) over one wavelength in permittivity 
e{z) (illustrated on top), for A=E— 6.5 cm = 3.16 nm and contrast Ase/ee = 0.1. The 
Casimir force profile [upper panel) is everywhere repulsive and ranges from zero for no 
effect at the permittivity wave maxima and minima to a maximum effect in the most 
slopping portions of the wave. Though most of the power is in the A/2 wave and 
bulk Fq modes, as expected from our analytic derivation in Section [3l a small residual 
remains, which is evident when the main Casimir wave modes are subtacted {middle 
panel). The residual mainly contains just the Fic and F^c modes. The Casimir force 
difference between normal and sign-reversed permittivity waves [bottom panel) closely 
follows the residual in the middle panel, showing twice its amplitude with matching 
relative Fic and F^c mode profiles. The main Fourier modes Fq and F2c do not change 
with the reversal of the permittivity wave phase and so vanish in the AF profile, however 
both of the other significant wave modes do fully. Thus the principle Fourier modes 
Fq and F2c are modulated with twice the frequency of the permittivity wave as we 
anticipated in Section [HI and the two residual modes Fic and F^c are modulated with 
the same frequency as the permittivity wave. 

All of the Casimir force profiles for small amplitude permittivity waves Asc/Sc ^ 0.1 
look qualitatively alike, being dominated by Casimir wave modes Fq and F2c. However 
the Casimir force profile at large permittivity contrasts does noticeable change its 
appearance as illustrated in Figure [9l showing an increased contribution from secondary 
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Figure 8. Casimir force profile F{z) for A=1E— 6.5 cni= 3.16 nm (upper panel), with 
principle wave modes Fq and i^2c subtracted (middle panel), and the force difference 
AF(z) between the normal and reversed permittivity waves (bottom panel), showing 
decomposition of profiles into wave modes Fic (dashed), F^^ (dot dash), and secondary 
residual F^c (long dash), for a permittivity wave of contrast Aeo/£o = 0.1 (illustrated 
on top), with ajc=2E16 rad/s and £0=6. 

Casimir wave modes Fi^ and F^^, and a small but here significant contribution from 
F^c {long dash). For larger amplitude permittivity waves, profiles of the residual force 
{middle panel), and force difference {lower panel) differ noticeably, indicating that the 
secondary wave mode F^c oscillates at twice the frequency of the permittivity wave with 
the principle wave modes Fq and F2c- 

For our numerical calculations the infinite integrals from (E]) are evaluated as finite 
angular integrals as for the multielement stacks in Section H] with sampling at 100 evenly 
spaced angular intervals. Integrals in ^ become negligible above the electronic vibration 
frequency and are cut off with a smooth function around lOOcJe- For general graded 
media, the integral in frequency vj is divergent, and so is smoothly truncated around 
the Compton frequency using a smooth Lorentzian, as we have described. 

Reflection coefficients at each ^ and vj were derived in two ways, using the recursion 
relation Eq. ( fT2l) starting many wavelengths from the domain of interest, and using the 
integral method from f|T6|) . For self-consistent results within a percent, each method 
requires many samples per unit of the dissipative length scale I/kq, 100 samples with 
the recursion relation (1T2|) . and 10 with the integral method ( TT6l) using Simpson's rule. 
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For the longest wavelength tested A=1E— 4 cm at the highest frequency ct7=lE20 rad/s 
for wavenumber = zu/c, 100ft:oA=3.3E7 recursion samples are needed or 10koA=3.2E6 
integration samples per wavelength. The integral method was iterated three times in 
the reflection coefficients R^s^_{z), and always gave good convergence showing changes 
no larger than 10% after the initial iteration. 

Both methods give nearly identical results, except that the recursion method is 
noisier, with noise comparable to signal for a relative permittivity fluctuation of about 
IE— 4. The integral solution gives consistent solutions even for relative permittivity 
fluctuations as small as IE— 8. The figures shown here were obtained using the integral 
method. 

The strength of the Casimir effect for a fluctuating static permittivity is 
approximated using f l25|) evaluating the integral with the permittivity model fl27|) for a 
nonconductor = and estimating the wave contrast by the derivative iic/io = Ae/e 
with Ae = {de/dec)Aec, which gives an analytic integral that is evaluated as 

For the parameters used in the Figure [HI the force has its most negative value for 
cos = —1 of F = — 1.31E10 dyn/cm^, a little larger than the most negative value for 
the numerical integrations from the figure. The small discrepancy arises because the 
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frequency cutoff in the numerical integrations is approximated by a smooth function 
rather than the sharp cutoff used in the analjd;ic integral from fl29|) . 

Figure fTOl shows the Casimir wave modes -Ffcc(A) {upper panel) and wave modes in 
the force difference AFfcc(A) {lower panel) as functions of wavelength A for a number of 
static-permittivity contrasts A6e/ee= E— 0.3, E— 1, E— 2, E— 3, E— 4, and E— 5. Each of 
the wave modes follows the power law predicted in Section |3l 
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Figure 10. Casimir wave modes for tlie force ^fcc(A) {upper panel) and the force 
difference AFfcc(A) {lower panel) as functions of wavelength A for different static- 
permittivity fluctuations Aec/^c showing the principle wave modes Fq and 
(everywhere coincident solid), residuals Fic and F^c (coincident dashed), and F^c {long 
dash). As wave modes overlap considerably, vectors are drawn to connect the successive 
Fq with Fic and Fic with F^c in each permittivity contrast group. Only Fic and F^c 
appear in the force difference. 



For the Casimir force {upper panel), the principle wave modes {solid) are 
separated by two orders 2 log]^Q(A£:/£:), so its amplitude goes like the permittivity 
contrast squared. The residual wave modes {dash) and F^^ {dot dash) always 
exhibit the same power and are exactly one order logio(A£:/£:) below the principle 
{solid) in each permittivity contrast group. The Casimir wave mode F^c {long dash) is 
also consistently reproduced, but everywhere two orders 21ogio(A£/e) below the main 
component in each permittivity contrast group. In every case, the bulk mode Fq follows 
exactly the wave mode F2c- In the difference spectra {lower panel) the principle Casimir 
wave mode and residual F^c disappear entirely. The difference spectra are composed 
of only the residuals Fic {dash) and F^c {dot dash), which exactly overlap and are 
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separated by two orders of permittivity contrast between groups. 

Comparing the vertical scales between Figures [6] to [10] shows that the Casimir 
force for a sinusoidal wave of small wavelength scale A=1E— 7 cm=l nm is almost a 
1000 times larger than what is found for the equivalent square permittivity wave or 
multilayer stack. The relative sinusoidal-wave strength grows considerably larger too 
at larger wavelength since the Casimir force for for sinusoidal waves goes only as 
whereas for square waves or multilayer stacks drops off as A^^ or A~^. 
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Figure 11. Casimir wave modes for the force Fkc{X) {upper panel) and the force 
difference AFkc{X) (lower panel) as functions of wavelength A for different plasma 
frequency fluctuations Awp/wp identifying wave modes as in Figure [TOl A plasma 
frequency of a;p=1.6E16 rad/s was used characteristic of the best normal conductors. 



Almost identical results are obtained with permittivity variations due to 
fluctuations in the electronic vibration frequency u^, so these are not shown. That 
solution is approximated from fl25|) with Ae = {de/duJc)A.uJc, giving 



F{z) 



TtHUcUc (gc - 11 

16A2c rpy2 



1)= 



Up 



Airz 

cos — 1 

A 



(30) 



Equating the force amplitude with that of (12 9 p with the same relative permittivity 
contrasts Aec/Ec = Auc/uc and solving numerically the transcendental equation shows 
that the forces are equal when ee=6. 038948 , very close to the static permittivity we 
use in our models. A slightly larger Casimir force is obtained with static permittivity 
fluctuations in materials of lower static permittivity, and with fluctuations in electronic 
vibration frequency in materials of higher static permittivity. 
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Casimir profiles for a fluctuating plasma frequency or Drude collision frequency 
resemble the profiles shown in Figures |8] and |9], except that the sign of the fluctuations 
is reversed for fluctuations in the Drude collision frequency. Figure [TT] shows the Casimir 
mode amplitude as a function of acoustic wavelength for a fluctuating plasma frequency 
with plasma frequency a;p=1.6E16 rad/s and Drude collision frequency 7p=lE14 rad/s 
for good normal conductors (with ee=Q and a;e=2E12 rad/s). The Casimir wave modes 
Fic and F^c (dashed) here actually slightly exceed the principle modes Fq and (solid) 
in the largest amplitude wave. The amplitude of the Casimir effect for plasma frequency 
fluctuations is proportional to the plasma frequency, and is little affected by changes in 
the other parameters. 

Similar numerical tests made using the Matsubara sum from ([3]) in place of the 
integral from ([6]), with reflection coefficients derived using the integral method (fT6l) . 
show no temperature effect in the Casimir force in permittivity waves, having tested 
with contrasts up to Ae/e=0.5 and out to wavelengths as long as A =E— 3.5 cm=3.16 
fim. We see no temperature effect whatsoever, perturbing any of the four parameters 
that enter the permittivity model (|26|) . 

6. Piezo ZP Energy Dynamics 

The practical utilization of zp energy requires that the Casimir-induced variations in 
a crystal couple effectively to an external electrical load. In this section we consider 
the piezo effect for the Casimir wave solutions from Eq. ( l25l) from Section [3] in both 
traveling and standing waves. The general energy transfer is described by assuming the 
conservation of an internal energy U (equivalent to the Gibbs free energy) in a stressed 
piezo crystal [32j 



using the summation rule over repeated indices in a term, with the usual stress Tij and 
strain Sij 3x3 symmetric tensors containing both longitudinal and sheared components, 
and the electric field Ei and electric displacement Di 3-element spatial vectors. Indices 
count coordinates i = x,y,z; Xi is the spatial coordinate (xi, 0:2, 0:3) = {x,y,z). 
For nonpyroelectric or magnetically affected piezo crystals and crystals that do not 
experience changes of state in the relevant ranges, the internal energy must be constant 
dU = 0, and the Casimir- forced internal stress couples exclusively to the crystal 
capacitance, so the total mechanical work T^dSij over one wave cycle must balance 
the capacitive energy gain —EidDi. 

Material parameters are taken to be linearly coupled connected by two constitutive 
relations [32]. We choose the most convenient pair for this work [33] 



with the constant impermittivity at constant stress /3^j a 3 x 3 real matrix, which is the 
inverse of the unperturbed permittivity matrix Sij for low frequencies where the electric 



dU = TijdSij + EidDi 



(31) 




(32) 
(33) 
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fields are important in the generation of capacitances tliat drive external currents. 
The elastic stiffness Cijke isa3x3x3x3 real tensor, which is the inverse of the 
elastic compliance tensor Sijke, CijpqSpgM = S(ij)6(^ki), for the Kronecker delta The 
piezoelectric constant hijk is a 3 x 3 x 3 real tensor, which enters symmetrically between 
the material response for the inverse piezo effect in (1321) and in the voltage response for 
the normal piezo effect in fl33|) . 

Stress-caused energetic processes do not give a capacitive transfer when the 
mechanical work done TijdSij over a cycle vanishes. In such cases the stress Tij and 
strain Sij vary in phase. Then the electric displacement Di and electric field Ei must 
also vary in phase for conservation of the total energy U from ( pTl) and for consistency 
with the constitutive relations fl32|) and f l33|) . Such zp energetics are conservative, so 
that whatever Casimir work is done in the compressive part of the cycle is returned in 
the elastic reaction, giving a Casimir force with no external effect. 

As described in the introduction, in a stack of alternating tunably conductive and 
piezo layers, the Casimir force might be equated to the material tension, which represents 
a static piezo response, relevant when the piezo layer is acoustically isolated from the 
rest of the stack and incommensurate in scale compared to the acoustic wavelength. 
The Casimir tension is on during the conductive or compressive part of the forcing cycle 
and off during the nonconductive elastic response, and thus maps out the area W as 
defined in Eq. ([2]) in the T-S phase space. As that calculation simply took the stress 
T to be proportional to the strain S and ignored the additional electrical displacement 
contribution D in ( l32l) representing the converse piezo effect, it needs to be increased 
by a multiplicative factor 1 + kp, equating the amplitudes of the total mechanical work 
and dielectric energies from (!3T1) . The piezo coupling constant kp is defined in Section 
6.31 of Mason 1966 [32j. It is typically small, kp ~ 0.1 in Mason's example. 

In general the Casimir force F{z,t) in a permittivity wave in a crystal acts to 
introduce coupled stress Tij, strain Sij, electric field Ei, and electric displacements 
Di, which depend upon the orientation of the natural crystal axes with respect to 
the acoustic wave z axis of the system. The attractive Casimir force per unit area 
is a traction force that acts internally in the crystal in the z direction and only as a 
function of z in our plane-parallel approximation, and thus must satisfy the force balance 
equation [3ll Section 2.B] 

d^u, dF{z,t) ^ 

with p the density of the medium, and Ui the particle displacement 3 vector. In 
addition the spatially varying component to must vanish due to Maxwell's equation 
in an acoustic wave for charge conservation, written for a nonconducting medium as 
V ■ {Di} = dDi/dxi = 0, denoting the vector of elements as {Di}, which reduces to 
its z component dD^/dz = when d/dx = d/dy = as we suppose everywhere in the 
interior of the crystal. Also Maxwell's vanishing curl of the electric field in an acoustic 
field V X {Ei} = 0, requires that the spatially varying transverse electric field in {E^j., Ey) 
also be zero. 
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Particle displacements define the strain tensor 

+ (35) 



2 \ dxj dxi 

which, for vanishing x and y derivatives everywhere in the interior of the crystal, gives 
the Casimir- force driving equations for the three nonvanishing strain elements S^z, S^z, 
and Syz, by taking the z derivative of ([3] 



^ dt^ dz^ dz^ ' ^ ' 

2piir = (37) 

for i = x,y. 

Using Eq. ( I33l) to eliminate the tension terms in the force equations ( l36l) and ( 1371) 
defines the strain for a given Casimir force F and electric displacement Di, for Szz 

P-^ - c?z.z^ = ^ (Fiz, t) - h^zM , (38) 
and for Sxz and Syz 

dt^ ^'^^^ dz^ ~ dz^ • 

The set of coupled linear differential equations are simplified by Fourier and 
eigenvalue analyses. Only strain waves and electric displacements at the frequency 
2v for the Casimir forcing are relevant, but because of the Casimir bulk mode in 
the driving force from ( l25l) . a general spectrum of wavenumbers kz are possible. 
We analytically continue by supposing the complex temporal variation in every term 
exp(?47rz/t) and Fourier transform with spatial z functions going over into wavenumber 
kz functions denoted by overbars. Then diagonalizing the matrix that multiplies the 
longitudinal/shear strain vector {S^^zikz), Syz{kz), Szz{kz)) from the left sides of Eqs. fl38|) 
and ( l39l) gives eigenvalue equations for at most three independent strain eigenvectors 

^(^^) = VI^2 {V^{kz) + hD{kz)) ; (40) 

'^zO '^z 

where the stain eigenstate is denoted S{kz). The Casimir forcing function is denoted 
F{kz) with ip a factor that projects the purely longitudinal z force into the strain 
eigenvector, and the projection of the shear electric displacement {D^, Dy) written D{kz) 
with its projection into the eigenvector contained in the single scalar coefficient h for 
the effective piezo force. The resonant wavenumber is defined fc^o = 47rz//cs for a sound 
speed Cs that contains the density p and the crystal stiffness from c^^^ for the eigenvector. 
The common temporal variation exp(i47rz/t) has been factored out of each of the terms, 
but possible temporal phase delays for the strain and electric displacement compared 
to the Casimir forcing may be contained as complex factors exp(z55') or exp(z5£i) in the 
functions S{kz) and D{kz). 

A transfer of energy between the mechanical and dielectric terms occurs if the 
integral of TijdSij (or equally —EijdDij) over a cycle is nonzero. As can be seen by 
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examining T-S phase diagrams, if the strain lags behind the stress, the integral TdS 
is positive and mechanical work is converted to dielectric energy. In some processes 
it is possible that the strain may lead the stress, corresponding to the transfer from 
an external electrical load into the quantum zp. The possible electric displacement and 
electric fields are defined by the boundary conditions, and by their coupling to the stress 
and strain from fl32|) and fl33|) . but with phase delays that depend upon the external 
circuit load. 

We examine the piezo energy transfer arising with the Casimir wave and bulk 
modes for a permittivity wave in a uniform medium from (l25l) of Section [31 A strain 
wave eigenstate may be driven as described in Eq. ( HOl) by the Casimir wave (standing 
or traveling) from fl25l) at wavenumber kz = |47r/A| if it is near resonance with the 
vanishing denominator at kz = kzo = 47rz//cs. Resonance thus requires that A = I'/cs, or 
that the sound speed for the acoustic wave is the same as that for the strain eigenstate, 
which has twice the frequency (and wavenumber). However high frequency acoustic 
waves with wavelength much smaller than the scale of the crystal do not normally 
couple capacitively anyway, as the spatially varying displacement field Di must be 
perpendicular to the wave direction and the electric field Ei parallel, and being always 
mutually perpendicular can do no work in the energy Eq. (13T1) . We can imagine certain 
exotic arrangements perhaps like what is described in the literature [35] that may allow 
piezo coupling to the main Casimir wave mode. However with such arrangements it may 
be difficult to avoid dissipating the acoustic wave too that drives the Casimir effect. 

The Casimir bulk mode for standing waves from (l25l) couples more straightfor- 
wardly. That mode is described as a spatial boxcar between uniform z-facing electrodes 
as in the drawing Figure [2l Charge buildup on the electrode surfaces defines a boxcar 
function too for the electric displacement Dz{t) giving sine functions in the kz Fourier 
domain for F{kz) and D{kz) in fHOj) . Thus the strain S{kz) is proportional the multiplier 
^2/(^20 ~ ^2) times a sine function in wavenumber space. Both functions are illustrated 
in Figure [T2l 




Figure 12. Wavenumber-space multiplier and forcing function for a strain response 
comparable in scale to the dimension of the crystal. 

If the resulting strain spatial function S{z) is small in magnitude compared to the 
electric displacement for the projected eigenstate \S{z)\ <^ 1-0(2) | over z, then the stress 
T{z) and electric field E{z) are determined solely by the electric displacement using 
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the constitutive relations and (ES]), and significant phase differences between the 
quantities are precluded, so no energy transfer possible. The total mechanical work 
done through the volume in a wave cycle is determined by the spatial-temporal integral 
of T{z,t)dS{z,t), which feels its determining contribution from the spatial-temporal 
integral of S{z, t) times the boxcar forcing. By the Fourier power theorem, the spatial 
integral can be equated to a wavenumber integral of the multiplier for S{kz,t) from 
fHOj) times a sine function squared. As long as strain variations are on the scale of the 
crystal, the integral over wavenumber of the product of the two functions is large and 
comparable to what would be obtained if the strain response were a pure boxcar. A 
significant strain response to longitudinal boxcar forcing on the scale of the crystal thus 
requires a small wavenumber k^o ~ corresponding to a wavelength for the eigenstate 
comparable to the size of the crystal. Characteristically phonon optical modes do exhibit 
zero wavenumber crossings corresponding to the Brillouin zone center at around 10 THz 
frequency in crystals [SU Chapter 3]. 

Perpendicular uniform surface electrodes say in x might be effective too, as electric 
displacement Dy{z) and electric fields Ey{z) constant in y may strongly couple to the 
induced strain field S{z). As the Casimir forcing is everywhere one directional, power 
generation must be relatively tolerant against incoherence in the standing waves due 
to partial or incomplete acoustic wave reflection, crystal imperfections, or acoustic 
attenuation effects. 

The voltage difference between electrodes is defined by the electric field, the order 
of which is determined by the piezo relation E ~ gT ~ —gF for g the piezoelectric 
voltage constant; characteristically g ~7E— 5 Volt cm/dyn for crystals. For a strength 
of Casimir bulk mode from Figure fTOl of around 5E8 dyn/cm'^, we estimate an electric 
field of £'=35000 V/cm, which gives a current of 28.6 Amperes per cm^ of electrode 
surface area for a nominal total power of 1 MW/cm'^, which is the product of the two. 

Though it has not been our main interest, stacks of alternating layers, the basis 
for our thought experiment, might give another practical design for a crystal power 
generator. A stack of alternating p- and n-doped piezo semiconductor layers are 
conductive throughout except in the depletion zones between the layers, which may 
be tuned conductively via electric or light fields. The design satisfies the need for 
fast and easy tuning as semiconductors operating differentially at p- to n- doped 
junctions exhibit a low excitation potential, and the Casimir force across the thin 
nonconductive depletion zone between conductive layers exhibits a strong piezo effect 
in many semiconductor materials. Fine tuning of the material parameters, doping 
concentrations, layer thicknesses, and frequency are needed to fix the thickness of the 
depletion zone. The light-travel distance c/z/= 30 fim for z/=10 THz is smaller than 
the size of multilayer stacks of interest, but larger than a useful layer thickness, but 
the stack need not be modulated in-phase. With slow in-phase cycling, an AC voltage 
across the total stack area and thickness is produced, whereas fast out-of-phase cycling 
feels differing contributions across the stack at any one time giving a more nearly DC 
voltage. 
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Figure 13. Available power per unit of crystal volume W as function of layer 
thickness X ^ £^ + £^ with £^ = i?_ in a stack with alternating static permittivity 
(a), electronic vibration frequency (b), and plasma frequency (c) for modulation at 
frequencies corresponding to the sound speed Cg = i^A=5E5 cm/s. Horizontal power 
levels (dashed) are for material compliance s=lE— 13 cm^/dyn. 

A stack of alternating layers cycled at a rate i' produces a Casimir forcing that 
cycles at twice that rate 2z/ as discussed at the end of Section HI Thus our estimate 
for the maximum power that might be extracted represented by the formula ([2]) needs 
u replaced by 2i/. Noting too the participation of the full volume in the process and 
supposing equal even-odd layer thicknesses = we obtain 

W = si^F^. (41) 

The formula ignores the small force difference AF ^ F between the two half cycles. 
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Figure 14. Available power per unit of crystal volume W as function of wavelength 
A in a static-permittivity wave of indicated contrast for modulation at frequencies v 
defined by the sound speed Cg — i^A=5E5 cm/s. Horizontal power levels (dashed) are 
for a material compliance s=lE— 13 cm^/dyn. 



Figure [13] gives estimates of the available power from ( HTll for modulating the 
static permittivity {panel (a)), the electronic frequency (b), and the plasma 
frequency (c), using the numerical computations from Section H] supposing the acoustic 
relation between frequency and wavelength z/ = Cg/A, and using a sound speed 
Cs=5E5 cm/s, characteristic of dense materials at least at frequencies below 2 THz, 
or A=2.5E— 7 cm=2.5 nm. The varying dispersion and possible phonon modes at small 
wavelengths may change greatly the appearance of the figure, but power levels should 
be representative. The available power goes like the Casimir force squared, and so 
drops quite rapidly with A, and also depends quite critically on the relative permittivity 
fluctuation going as the 4th power. 

Estimates for the plasma frequency modulation in (c) are similar in their maximum 
values, but form a relatively tight group for the range of possible plasma frequencies. For 
an ideahzed perfect conductor with plasma frequency above a;p=lE18 rad/s, the curves 
show little change, and represent an available power of about 1 MW/cm^ with £o=5 
nm (A=10 nm=lE— 6 cm), for the conservative piezo compliance s=lE— 13 cm^/dyn, 
consistent with the estimate given in the introduction. However power levels drop off 
rapidly when the plasma frequency contrast is less than 1E18 rad/s as with normally 
conductive materials, which exhibit plasma frequencies only as high as 1.6E16 rad/s. 

Figure [T3] shows power per unit volume W from (l^Tjl for the main Casimir modes, 
the bulk mode Fq or main wave mode as they are of equal amplitude, as a function 
of the acoustic wave wavelength for static-permittivity fluctuations of differing relative 
strength. The figure illustrates the high power levels that are possible even with a modest 
relative permittivity wave amplitude of A£/e=lE— 3 and shows the much greater range 
into longer wavelengths for the permittivity waves compared to multilayer stacks. Such 
amplitudes are easily below normal crystal damage thresholds [20] . 
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8. Mode Growth and Attenuation 

Casimir modes should arise in any material, not just a piezo crystal, and may feed back 
producing either growth or attenuation in the principle acoustic wave. An acoustic wave 
exhibits an associated permittivity variation due to the acousto-optic effect, which is 
determined phenomenologically by the Pockels second-order nonlinear response due to 
the piezo-induced electric field caused by an acoustic wave of pressure amplitude P 

^ = epsP, (42) 

with s the material compliance relevant for the propagation direction and shear. The 
photoelastic coefficient p is a dimensionless parameter for the material in the wave 
direction for an unsheared or sheared stress for longitudinal or transverse waves, 
respectively; characteristically p ~ 0.1 and usually positive. 

The Casimir mode Fic must always couple linearly to the principle acoustic wave 
as it is of same frequency and wavenumber. For materials of positive photoelastic 
coefficient p, the acousto-optic effect from fH2|) exhibits a pressure variation correlated 
with the induced permittivity wave b{z), giving a Casimir wave mode Fic with the 
same sign as the pressure in the acoustic wave and so amplifying. But in materials of 
negative photoelastic coefficient p < 0, such as GaAs or Ruby, the Casimir wave mode 
Fic is dissipative. In each acoustic wave cycle for standing or traveling waves at a given 
spatial location, the Casimir mode Fic is generated from quantum zp energy and forces 
the principle acoustic wave. All of the energy in the Casimir mode Fic feeding into the 
pressure wave amplitude P gives the maximum growth rate 

Tic = = 27rz/-^. (43) 

Taking Fic from f l29|) with P from fH2|) . and using the static-permittivity fluctuation 
for the permittivity variation e — )■ £e with the result from Figures IHlfTOl that Fid Fq — 
A^e/e^e, we obtain 

Tr'^hWf.uJcpse^J'^v^ ( IS.ec 



Tic 



32c ci V e, 



0.216 GHz 



3 /Aec 



(44) 



now using the acoustic wave frequency v rather than its wavelength A in the formula, 
with sound speed Cg = Az/, and using characteristic parameters for the rightmost 
equality p=0.1, s=lE— 13 cm^/dyn, Cs=5E5 cm/s, a;e=2E16 rad/s, a;c=7.76E20 rad/s, 
the Compton frequency, and £:e=6. The growth rate represents an amplitude e-folding 
in a distance Cg/F for the often-used attenuation factor a = — (201oge)F/cs dB, giving 



in units of decibels per cm. For a wave of frequency u = 1 THz and high amplitude 
A^e — £e, we obtain a maximum growth rate of 28 dB/cm. The growth rate drops off 
rapidly with decreasing acoustic wave amplitude going as (A^e/^e)^- 
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Unless it is an phonon eigenstate, the primary Casimir mode will be 
strongly dissipated by interference from out-of-phase driving by modes in adjacent 
wavelengths of the principle acoustic wave, leading to a rapid exponential attenuation 
rate. Though of frequency 2i>, the mode is coherent with the principle acoustic wave 
and might exchange energy effectively with it if there is a sufficient nonlinear interaction 
between the two. Coherent acoustic modes may exhibit nonlinear interactions similar to 
coherent electromagnetic states, which exchange energy easily between overtones, and 
generally satisfy a principle of energy equipartition [36]. If there is a strong coupling, 
the attenuation rate for the principle acoustic wave may be as much as half the 
mode amplitude per acoustic wave cycle to give an attenuation rate or negative growth 
rate in the principle acoustic wave a factor (Aee/f^e) larger than the growth rate for the 
smaller Fic mode from ( 143|) - ( l45l) . Such attenuation would go like and be nonlinear 
as it depends upon the mode amplitude. 

If the Casimir mode F2C is a longitudinal phonon eigenstate, it may similarly 
affect growth in the principle acoustic wave of comparable magnitude. Thus it may 
be desirable to match frequencies and crystal orientations in materials where both of 
the main Casimir modes with frequency 2i> are phonon eigenstates: the bulk mode with 
wavenumber /c^ ~ and wave mode of wavenumber twice that of the acoustic wave. 
Further discussion on the growth or attenuation of the principle acoustic wave that 
might be caused by Casimir modes is material dependent and discussion of possible 
nonlinear interactions seems beyond the scope of this paper. 

Noncrystalline substances do exhibit a limiting high-frequency attenuation that 
increases faster than z/^ and appears to be close to a i/^ form [371 Figure 10]. Though such 
measured THz attenuation factors are much larger than ours, these are for noncrystalline 
substances, which are known to exhibit a much stronger acoustic attenuation than what 
is seen in pure crystals [381 ESI HQ]. Attenuation effects exhibit a strong temperature 
dependence below 100 GHz, but become temperature independent above that frequency, 
as we expect for Casimir modes. Measured attenuation effects in noncrystalline materials 
are found to have a nonlinear amplitude dependence at least at 10 GHz frequency |39j . 

9. Conclusion 

We have studied the Casimir effect in multilayer stacks and in graded-permittivity media 
both analytically and with numerical models. In a stack of alternating high- and low- 
permittivity layers, an attractive Casimir force is produced throughout the volume, 
which follows closely solutions for a three-layer sandwich, with thin- and thick- central- 
layer £0 limiting power-law forms that go like l/^g and l/i^, respectively. The Casimir 
force switches power laws at fairly small scales depending upon the permittivity contrast, 
as at thickness ~ 7 nm for usual parameters. The compressive force is only slightly 
smaller across the high-permittivity layers than across the low-permittivity layers. By 
modulating the permittivity in the layers, it is possible to cycle the Casimir force on and 
off through the volume and harvest energy in each cycle by the piezo effect. The concept 
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might be implemented in a stack of alternating p- and n- doped piezo semiconductor 
layers, which are conductive and might rapidly tune the conducting state of depletion 
zones between the layers by an electric field or light. 

Similarly a short wavelength permittivity wave in a uniform material produces 
large Casimir half-wavelength and bulk expansive modes. The Casimir effect in a 1 nm 
acoustic wave is about 1000 times stronger than that found in a stack of alternating 
layers of the same permittivity contrast and scale of variation, and extends to much 
larger scales decreasing only as The Casimir effect in a wave does not feel a 

natural frequency cutoff as arises across a uniform layer due to retardation extinction, 
but feels contributions out to the frequency cutoff for the zp interactions, which we take 
to be the Compton frequency. Smaller amplitude Casimir modes arise too as higher 
wavenumber overtones of the principle wave with amplitudes proportional to increasing 
powers of the relative permittivity contrast. Unlike for a multilayer, the Casimir effect 
in a permittivity wave exhibits no temperature dependence at least within the scope of 
our numerical models. 

Casimir modes induced by THz acoustic waves in a piezo crystal may couple 
effectively to an external electrical circuit, in effect taking energy out of, or returning it 
to the quantum zp, as we have shown through the piezo thermodynamic formalism. A 
significant energy transfer is possible for standing acoustic waves of frequency u through 
the Casimir bulk mode, but only if matched to a low wavenumber eigenstate at frequency 
2i/, like a phonon optical mode. Acoustic waves might be introduced and coherently 
maintained in different ways, like by shining monochromatic microwave light on one or 
opposite piezo crystal surfaces through transparent electrodes [H]. The standing wave 
coherence that is possible for nm waves might be greatly limited by surface or crystal 
imperfections, or other sources of acoustic attenuation. However since the Casimir bulk 
mode is expansive, the voltage effect is single signed and so should add constructively 
through the volume with partially coherent standing waves for a partial ac or dc voltage 
difference across the crystal. 

We estimate that significant power may be derived of maximum order a MW/cm^ 
for 1 nm=lE— 7 cm wave even with a modest relative permittivity contrast of IE— 3. 
A significant effect persists even out to 100 nm acoustic wavelengths due to the slow 
falloff of the Casimir force with wavelength. 

Casimir modes may be important in the attenuation or driving of hypersonic (10 
GHz - 10 THz) acoustic waves through nonlinear or direct linear feedback, as we have 
discussed. 

Author's Note 

It is the author's expressed intent that the guidance, concepts, and designs described 
remain unpatentable but in the public domain. It is his belief that energy derived from 
the quantum zp is God given for the free and beneficial use by all on His good earth. 
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Appendix A. Limiting Solutions for Multilayer Stacks 

Power laws arise for the Casimir force in the thin- and thick-layer limits for certain forms 
of permittivity variations. In this appendix, we parallel the Lifshitz derivation [121 EI] 
for the Casimir force integral from ([6]) for permittivity variations defined by contrasts 
in each of the parameters contained in ( 1271) . 

A rearrangement of the Casimir force integral in facilitates evaluation of the 
force in the middle layer. Using the new variable x = 2/«o^o = (£o(0'C^ + ^^)^^^^o/c to 
eliminate kq under the integral, and substituting for the integration variable w using 
dx^ = {Ml/c^)dw'^, gives 



denoting layers by subscript number with Eq^ and kq for the middle layer. The force 
integral in (lA.ll) is dominated by the single value x — 1 due to the tradeoff between 
the weighting factor under the integral, which increases with x, and the decreasing 
exponential exp(— 2ko£o) = exp(— x), which determines the reflection terms Rts\ in 



through the products of Fresnel reflection coeflicients from Eq. (iTTl) . 

As Lifshitz points out [12j, when the middle layer is sufliciently thin £o ^ c/cUg, the 
dominant wavenumber kq ~ l/^o is higher than the wavenumber for electronic vibrations 
^ ojc/c. The permittivity between layers changes little at frequencies uj^ 
as is evident in Eq. ( !27l) . so the wavenumber between layers Kj changes little too. 
The jump coeflicients rs± from ([8]) contribute little to the integral, and s polarization 
effects and Rs can be ignored. Also dependencies upon n in the p polarization jump 
coefficients rp± from (jO]) or in Rp divide out. Taking contributions from only the nearest 
jump in the Fresnel reffection terms, and dropping the product of reffection coefficients 
from the denominator of Hp in (j4]), as it introduces a fourth-order effect in the relative 
permittivity change compared to the main second-order effect, we obtain the Lifshitz 
thin-layer Casimir force formula written for permittivities symmetric about a middle 



using a compact notation to denote the permittivity contrast = {eq — e±i)/2 = 
(£__ — £+)/2, and the average e = {eq + e±i)/2 = {e^ + e+)/2 for layers numbered 0, 
±1, even — , or odd +. As the permittivity e(^) does not depend upon x, the x integral 
factors out in the derivation as a Gamma function x^^~^dx = 2. 

Supposing that the permittivity contrast arises due to static-permittivity 
ffuctuations Ae = {dE/dEe)AEe for permittivities deflned in (1271) . we obtain an analytic 
integral in (]A.2p with the solution 




(Al) 



J = (Eq. (4.18) m) 




(A2) 




(A3) 
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assuming a nonconductor oOp = for simplicity. The thin-layer limiting force with a 
fluctuating electronic vibration frequency As = {d6/duJc)Aue similarly gives 

Corresponding limiting behavior with varying plasma frequency or Drude parameter is 
less obvious, as suggested in the numerical solutions described in Section HI 

When the middle layer is sufficiently thick, the main contribution to the Casimir 
force in flA.ip comes from low frequencies, and wavenumber changes from layer to 



layer cannot be ignored. If static-permittivity variations exist, then it seems a good 
approximation to suppose that the permittivity is constant in frequency and well 
approximated by its static value = 0) . A critical integral substitution allows 
all of the variables of the problem to factor out and renders the integral essentially 
constant over parameter ranges of interest. Such a substitution is suggested by the 
behavior of the jump coefficients Tg- in the low-frequency limit, which are equal to the 
wavenumber contrast between layers r^- = Ak/k. Differentiating k with respect to e in 
the wavenumber definition in ([5]) allows us to translate the wavenumber contrast into a 
static-permittivity contrast as 

Ak_[ ^y^y AsiO) ^ fi^y AsiO) 



2cK j e{0) \xj e{0) ' 

where x = 2ko£o — 2k£o for small fluctuations Ak/k, and deflning a new variable 
= ^^o-\/^(0)/c. By substituting the new variable v for ^, the wavenumber contrast 
becomes consolidated into a permittivity contrast with dimensionless coefficients defined 
by the variables of integration. Such a substitution indeed factors out all spatial scale 
dependencies and simplifies the integral in flA.ip . giving 

' Vm4 V £(0) J ^ ' 

with the coefficient fi approximately constant over a wide range of static permittivity 
contrasts. Lifshitz has shown this simplification to an equivalent expression for 
noninfinitesimal permittivity contrasts, which leads to a nonanalytic exponential-type 
integral he denotes 0dd(£i(O)/£o(O)), which relates to our coefficient fi = (7r^/24O)0dd(l + 
2A£:(0)/£:(0)), writing the argument of 0dd for small static-permittivity fiuctuations 
A£:(0) ^ £^(0). He evaluates the integral numerically and obtains (p^d = 0.35 for 
ei(0)/£o(0) > 0.3 (Eq. (4.22) [21]), which gives fi = 0.0144 for relative static- 
permittivity fiuctuations Ae{0)/e{0) > —0.35. From the Lifshitz more general 
formulation, when the sandwiching layers are nonconductors with large permittivity 
contrast {Aed ^ £c or pure conductors, the Casimir force is approximated as in ([1]) 
discussed in the introduction. 

Only as long as the main permittivity fiuctuations have a low- frequency component, 
as in fiuctuations of the static permittivity ^e, the plasma frequency cOc, or the Drude 
collision frequency Fp does the limiting thick-layer solution form (1A.6I) apply. It does not 
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apply with permittivity fluctuations produced by fluctuations of the electronic vibration 
frequency as illustrated in the numerical examples from Section |H 

The cross-over thickness between the two power laws iQ is defined as P{{^{)) = 
F){iQ). Using (lA.Sp and (1A.6P for static-permittivity fluctuations, we obtain 

= 1.448^, (A.7) 

supposing that a permittivity contrast arises due to static-permittivity fluctuations with 
A£(0) = Aec and e(0) = ^e- For the case in Figure [5] with a;e=2E16 rad/s and ee=6, we 
obtain £^)=3.62E— 7 cm or A()=7.24E— 7 cm=7.24 nm=E— 6.14 cm, in good agreement 
with the intersection of lines seen in the figure. The cross-over distance decreases 
with increasing static permittivity, suggesting a smaller cross-over distance in more 
conductive media. 



Appendix B. Computational Considerations with a Graded Medium 

Calculations of the Casimir effect in a medium with a small amplitude graded ID 
permittivity variation are based upon the force integral Eq. ([6]) with (jlj), with Fresnel 
reflection coefficients defined by the integral ( fT6i) with (HM and ( |T5l) . For numerical 
calculations the spatial integral in f|T6|) is best defined in the specific way 

R^s^yiz) = T^.)_(^) + e-2/^>(^')'^^'C^.^_, (B.l) 

where 



Tfs^_{z) = I e^^l W^")-'^(^))'^^"r^.^_(/) ( 1 - B!^.y{z'f ) d/, (B.2) 
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now including lower bounds on the integrals z^, for the integrating factor and z\y for the 
main integral. The lower bounds are arbitrary, and ^i^"^- is a constant of integration. 
We have taken the external exponential factor under the integral, which is a form more 
suitable for numerical calculations, as the wavenumber difference helps avoid very large 
exponential factors at the limits of normal computer precision that arise in integrals of 
the wavenumber Integrals over a full long wave A=l /im are of order 1E6 at the 

highest frequencies Ct7=6E20 rad/s, which would lead to cancelling exponential factors 
of order exp(lE6). 

Periodic solutions of wavelength A require the condition i?^s^_(zc) = -R^s^_(zc + 't-A) 
for any origin z^ and integer n. The condition uniquely determines the constant of 
integration since the solutions must be periodic, as we have described in Section 

[3 Substituting with f lRT]) then gives, with careful handling of the limits of integration 



in the sums and differences of integrals that arise, the expression 

Using this form to represent the constant of integration avoids other possible occurrences 
of large exponential factors. The term is the offset in fIB.ip . which may be directly added 



(B.3) 
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to the integrals T^s-j ^-z) to give periodic solutions for the Fresnel reflection coefficients 
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